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O , 

. We have proposed a general method for finding the exact analytical solution for the 

CN . 

H ' multi-channel curve crossing problem in presence of delta function couplings. We 

have analyzed the case where a potential energy curve couples to a continuum (in 
energy) of potential energy curves. 
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Ph: I. INTRODUCTION: 

^ ■ Nonadiabatic transition due to potential curve crossing is one of the most important mecha- 



nisms to effectively induce electronic transitions jl| . This is a very interdisciplinary concept 
and appears in almost all fields of physics and chemistry, even in biology [2]. The com- 
mon examples are variety of atomic and molecular processes such as atomic and molecular 
^ I collisions, chemical reactions and molecular spectroscopic processes. Attempts have been 
' made to classify organic reactions by curve crossing diagrams . There are also examples in 

o : n 1 

0\ ', nuclear physics [4|, l5| . Two state curve crossing can be classified into the following two cases 
O ■ 

according to the crossing scheme: (1) Landau-Zener (L.Z) case, in which the two diabatic 
^ . potential curves has the same signs for the slopes and (2) non-adiabatic tunnelling (N.T' 
\ case, in which the diabatic curves has the opposite sign of slopes. In our earlier paper p] 
paper we have considered two diabatic curves, crossing each other and there is a coupling 
between the two curves, which causes transitions from one curve to another. This transition 
would occur in the vicinity of the crossing point. In particular, it will occur in a narrow 
range of x, given by 

Vi{x)-V2{x)c:^Vu{Xc) (1) 

where x denotes the nuclear coordinate and Xc is the crossing point. Vi and V2 are determined 
by the shape of the diabatic curves and Vu represent the coupling between them. Therefore 
it is interesting to analyse a model, where coupling is localized in space near Xc- Thus we 
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Configurational Coordinate 



FIG. 1: The multi-channel curve crossing problem. 

put 

V^2{x) = K^5{x-x,) (2) 

here Kq is a constant. This model has the advantage that it can be exactly solved [gI. 
In this paper we apply this model in the case of general multi-channel problems (Fig. [1]), 
since this is very important in view of the applications to various practical problems in 
chemical dynamics beyond simple diatomic systems. Using model potential in which one 
inear diabatic potential intersects many constant diabatic potentials, Demkov and Osherov 
~\ showed that the state to state transition probabilities can be exactly calculated by a 
product of the Landau-Zener probabilities at each crossing. A more general model, in which 
one set of parallel linear diabatic potentials intersect with another set of parallel ones, 
has been analyzed by Child ^. Korsch jo^ reported multi-state curve crossing problems 
by the Magnus approximation, and later Nakamura [l| analyzed the problems by complex 
phase-integral technique. In the following we give a general procedure for finding the exact 
quantum mechanical solution for multi-channel curve crossing problem by extending our two 
state model 
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II. TWO STATE PROBLEM 



The time dependent Schro dinger equation (we take ^ = 1 here and in subsequent calcula- 
tions) 

(3) 



i — = H^{x,t). 



dt 



H is defined by 



where HJx) is 



and 



H 



HAx) 



Hi{x) V{x) 
V{x) H2{x 

1 



2m dx"^ 



+ VAx). 



(4) 



(5) 



(6) 



ipi{x,t) and 'i{j2{x,t) are the probability amplitude for the two states. We find it convenient 
to define the half Fourier Transform "^(u) of \E'(t) by 



Jo 



(7) 



Half Fourier transformation of Eq. ([3]) leads to 



u-Hi V 
V U-H2 



^i(O) 

^^2(0) 



(8) 



This may be written as 



^{uo) =iG{uj)^{Q). (9) 
G{(jj) is defined by {u — H) G{uj) = I. In the position representation, the above equation 



may be written as 



where G{x,xq;u) is 



'^{x,u)=i / G{x,xo;u)'^{xQ,u)dxo, 
J —00 

G{x, xq] io) = {x\{uj — H)^^\xo). 



Writing 



G{x, Xq] uj) 



G^^{x,xo;uj) G^2i^,xo;uj) 



^ G2iix,xo;uj) G22ix,xo;uj) 



(10) 



(11) 



(12) 



and we can write 



G^^{x,xq-uo) = {x\[uj- Hi-V{uj- H2)-^V]-'\xo). 



(13) 



This expression simplify considerably if y is a delta function located at Xc- In that case V 
may be written asV = KqS = Kq \xc) {xc\. Then 



G^^{x,xq]uj) = {x\[uj -Hi - KqGI(xc,Xc;uj)S] ^Ixq), 



where 



(14) 



(15) 



Gl{x,Xo;u) = {x\{uj - H^) Vo), 
Now we use the operator identity 

{uj-Hi-K^,Glixc, Xc; iu)S)'' = {uj-HiY^+{uj-H{)-^KIgI{xc, x,; uj)S[uj-Hi-KIgI{x,, x,; uj)SY\ 

Inserting the resolution of identity / = j'^^dy \y) {y\ in the second term of the above 
equation, we arrive at 



G]^^{x, Xq] Uj) = G^{x, Xq] U>) + KqG^{x, Xc] Uj) X G^l^^c, ^c] ^^)Gll{Xc, Xq] u) . 

Putting a; = Xc in Eq. (fT6l) and solving for Gii{xc, Xq; u>) gives 

Gu{xc,xo;u}) = G°{xc,xo;u}) x [1 - KqGi{xc,Xc;u)GI{xc,Xc;uj)]"^. 
This when substituted back into Eq. f|T6|) gives 



G'ii(x, xo; uj) = Gl{x, xq; uj) + 



KqGI^X, Xc] Uj)G'^{Xc, Xc] Uj)G1{Xc, Xq] Uj) 



1 - K^Gl{xc, Xc, uj)G^{xc, Xc, uj) 
Using the same procedure one can get 

'^12[-^,-^0,^) - l-K-^G'>ixc,Xc;oj)G\>{x,,x.;uj)- 

Similarly one can derive expressions for ^22(^)^0;^) and G2i(^>^o;^)- 



(16) 



(17) 



(18) 



(19) 



III. FORMULATION OF MULTI-CHANNEL CURVE CROSSING THEORY 

We start with the time dependent Schrodinger equation for a three state system (we put 
h = l) 



. d 

''Ft 



i/ji{x,t) 



\ 



i^2{x,t) 
\i)^{x,t) j 



( 



Hnix) (x) Vnix] 



\ 



Voa(x) Hooix 







iji{x,t) 

1p2{x,t) 

\ V^i{x) H^z{x) ) \ ^lJz{x,t) ) 



(20) 
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The above equation can be splitted into 



and 



^dij2^ = V2,ix)M^,t) + H22{x)ij2i^,t) (22) 



and also 



^d^h^ = Vs,{x)Mx,t) + H,,{x)Mx,t) (23) 

We define half Fourier transform \E'(a;) = '^{t)e^'^^dt and Fourier transform of the above 
three equations gives 



- z?/'i(x,0) +u;^i(x,cj) = Hii{x)il)^{x,uj) + Vi2{x)il)2{x,u) + Viz{x)il)^{x,u) (24) 

and 

u^p^ix.uj) = V2iix)ij-^{x,uj) + H22ix)ip2ix,uj) (25) 

as ip2{x, 0) = 0. Also 

u'(p^{x,u) = V3i{x)ij^{x,u) + H-i3{x)^^{x,u) (26) 
As ?/'2(a^,0) = 0. We will substitute -02(2^)^) S'^d '03(x,a;) in Eq. ([2] 



^2(^,0;) = [u; - //22(x)]-V2i(x)^i(x) (27) 

and 

^X.U) = [U- i/33(x)]-V3i(x)^i(x) (28) 

Eq. (1241) becomes 



[a; - iJn(2:) - Fi2(x)[u; - /f22(x)]^ ^21(0;) - V,3{x)[u - if33(x)]- V3i(x)]V^i(x, a.) = i^i{x, 0) 

(29) 
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The above equation simplify considerably if the couplings are assumed to be Dirac Delta 
function, which in operator notation may be written as V21 = V12 — K2S — -ftr2| 2^12) (2^12! 
and V31 = V13 = KzS — A'3|xi3)(a;i3|. The above equation may be written as 

[{oj-Hu{x))-Kl5{x-xi2)Gl{xi2,xi2]i^)-Kl5{x-xi^)Gl{xi^,xir^]uj)]i}^{^ = #i(a;,0) 

(30) 

or 

[{io - Huix)) - KlS{x - xi3)G2{xi3, X13; u)]iJi{x, u) = i'ilJi{x,0) (31) 

where, 

Hi2{x) = Hu + Kl5{x - xi2)G\{xi2,xi2]uj) (32) 
The solution of this equation may be expressed in terms of the Green's function G'(x, Xq\ u) 

[{uj - Hi2{x)) - Kl5{x - xi^)Gl{xi^, xi^] uj)]G'{x, xq; uj) = iS{x - xq) (33) 

Also 

00 

'ipi{x,u) — i j dxoG'{x^xo;u;)ipi{x,0) (34) 
—00 

Here G'{x, Xq] cu) is 

G'{x,xo;u;) ^ {x\{u} - Hi^y^lxo) (35) 

where 

His = Hu{x) + Kl5{x - xi3)G^(a;i3,xi3;u;) (36) 
In the previous section we have calculated Gii{x^xq]uj), which is redefined below 

Gii{x,Xq]uj) = {x\{lli - H^y-^lxo) (37) 

The Guix,XQ;u!) can be written 

r (r r ■::\-r^(r r ■ : :\ I KlGl{x,Xi2,i^)Gl{xi2,Xi2;Uj)G\{xi2,X^-u) 

Gn{x, xo, - G,ix, xo, + , _ k2gU,^^,xi2;u;)G%xi2,xu;u;) ^^^^ 
In this case, one can find the Green's function Gii{x,Xo;uj), using the method as we have 
used before. 

r'(r r -.A-r (r r ■ : A I ^iGajx, Xn] u;)Gl{xi3, X13; u;)G,,{xis, Xq; Uj) 

Lf [X, Xo, UJ — Lrii[X, Xq, LO H ^ r [6\) 

1 - KiGu{xi3, xi3;u;)G'i{xi3, X13; u;) 
So in this way one can add the effect of 4-th state, 5-th state etc. one by one and calculate 

Green's function for A/^-channel problem. 



A. Mult i- channel curve crossing problem: an analytically solvable continuum 

model 

Here we start with A^-channel generalization of Eq. (IMl [251 1261) 

N 

- z^i(a;,0) + t^^i(x,cj) = Hu{x)'i/j^{x,uj) + ^ Vie{x)'i{j^{x,uj) (40) 

e=2 

and 

u'iIj^{x,uj) = Vsi{x)'il)^{x,uj) + Hss{x)i!^{x,uj) (41) 
From the Eq. (14T1) we calculate '?/'g(x,a;) 

^,(x,u;) = [cu - if,,(x)]-V,i(x)^i(x,u;) (42) 
So Eq. fj40|) can be written as 

N 

- #i(a;,0) + ujt/j^{x,uj) = Hii{x)'iIj^{x,u) + ^ Vie{x)[u - Hee{x)]''^Vsi{x)'ilj^{x,u}) (43) 
we take Ki(a;) = Vis{x) = ks{x)6{x — x^). The above expression becomes 

_ _ Af _ 

— itpi{x, 0) + uJ^|J-^^{x,u) = Huix)ip-i^{x,uj) + ^ A;£(x)^5(x — a:£)G°(xe, x^; u)iJj-^^(x,uj) (44) 

£=2 

where 

G°(Xe,Xe;u;) = {Xe\[uJ - Heeix)]-'^\Xe) 
oo 

= E (a;£|</>n)[^ - //£e(a;)]"^(0nk£) (45) 

n=0 ^ ' 

oo 

= E (^^(X,)[.;-/7,,(X)]-Vr(^e) 
n=0 

In the above expression, e takes discrete values. Now 

TV 

[{u - Hii{x)) - ke{xY5{x - Xe)G^^{x^,x^]u)]^jj^{x,u) = i'ijji{x,0) (46) 

E=2 

The solution of this equation may be expressed in terms of the Green's function G"{x, Xq] u) 

N 

[{U - Hii{x)) - ke{xY5{x - Xe)Gl{Xe, Xe] Uj)]G"{x, Xq] u) = i6{x - Xq) (47) 

£ = 2 
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Here G"{x, Xo;u>) is 



N 



G"{x,Xo;Uj) = {x\[{uj - ifii(x)) - E ke{xY5{x - Xe)G'^{Xe,Xe]Uj)] ^|xo). AlsO 



e=2 



'ijj-^{x,uj) = i / dxQG"{x,Xo;u)'ijji{x,0) (48) 



In the continuum limit (here e changes continuously), the expression for G"{x, Xq; u) can be 
written 



G"{x,xo;uj) = {x\[{uj — Hii{x)) — J dek{x,ey6[x — x{e)]G^{x{e),x{e);uj,e}] ^\xo) (49) 



In a very special case where all states couples at one point to the first one i.e. x{e) = a 

oo 

G"{x,Xo;uj) = {x\[{uj — Hii{x)) — 6{x — a) / dek{x,eyG^{a^a;u,e)]~^\xo) (50) 



The above expression is very general and applicable to any potential. We can write Eq. fISU]) 
in a simplified form 

G"{x,xo;u) = {x\[{uj - Hn{x)) - 6{x - a)V{x,u]-^\xo) (51) 

where 

oo 

V{x,uj) = J dek{x,eYG^[a,a;u,e] (52) 







Special case 



In the following, we discuss the case where one potential (general) is coupled to continuum 
(in energy) of constant potentials at a point. Green's function for uncoupled potential is 
given by 



1 / TTl 

In the above expression e varies continuously from to oo. 

Analytical form of V{x,E) depends on the functional form of k{x,e). For our calculation 
we take 
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kix, = ad ^^^"^ 'K -it (54) 
V m 

So 

oo 

V{x,uj) = J dee~T^ = ahu (55) 



The Eq. flSTl) now becomes 

G"{x, Xq; uj) = {x\[{uj — Hii(x)) — S{x — a)ahuj]~^\x(i) (56) 
The above expression can be evaluated exactly. 



IV. CONCLUSIONS: 



We have extended our two state model to deal with general one dimensional multi-channel 
curve crossing problems in presence of delta function couplings. We have analyzed the case 
where a potential energy curve couples to a continuum (in energy) of potential energy curves. 
Nonadiabatic tunneling in an ideal one dimensional semi-infinite periodic potential system 
is also analyzed using our model. The same procedure is also applicable to the case where 
iS is a non-local operator, and may be represented by 5* = |/) Kq {g\ , f and g are arbitrary 
acceptable functions. Choosing both of them to be Gaussian will be an improvement over 
the delta function coupling model. S can also be a linear combination of such operators. 
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